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Abstract 

New classes of classically integrable models in the cosmological theories with a scalar 
field are obtained by using freedoms of defining time and fields. In particular, some 
models with the sum of exponential potentials in the flat spatial metric are shown to be 
integrable. The model with the Sine-Gordon potential can be solved in terms of analytic 
continuation of the non-periodic Toda field theory. 



After the original inflational scenario was presented, various models have been pro- 
posed to obtain a natural inflation such as the chaotic inflation and the modified Brans- 
Diche theories. In all of these scenarios, scalar fields play an important role. In other 
words, in the early universe, the dynamics seems to be governed by the gravity theory 
coupled with a scalar field. 

As for classical solutions, it is known that we can easily get them once we start by 
assuming the time variation of the scale factor a: The classical solution of a scalar field 
and the potential are expressed as functions of a and its derivatives so that in principle 
the potential is expressed as a function of a scalar field[l]. Moreover, equations for gauge 
invariant variables for scalar perturbations can be reduced to a second order homogenious 
differential equation whose coeficients are solely determined by a and its derivatives. By 
solving this eauation, several analytic solutions for the inflation models are derived[l]. 

In this paper, we present a method to get classically integrable models for gravitational 
theories coupled with a scalar field in the case where the dilatational mode and scalar 
field do not have spatial dependence. A method to get integrable models was proposed in 
Ref. [2], but unfortunately the integrable models were very few. Here, we show that with 
a slight modification of their approach, we can get integrable models for the larger class 
of models. In particular, it will be shown that models with the exponential potential 
corresponding to the Brans-Dicke theory are classically integrable. We also find the 
Sine-Gordon theory coupled to gravity can be solved in terms of Toda theories. 

The action of the gravity coupled to a scalar field is given by 

S = J d'x^i-J-^fi^R + ^{TWgTdfidv* - [/($))], (1) 

where /(<&), T($) and U(&) are functions of <3>. The standard choice of the action 
is / = 1 — £$ 2 ,T = 1 (the choice £ = | is called minimal coupling). The modified 
Brans-Dicke theories correspond to the choice / = $, T = cj($)/$ 2 , where the ordinary 
Brans-Dicke theory is given by to — 1, and other choice is called the extended Brans-Dicke 
theories[3]. For T($) = 0, the scalar field is just the auxirialy field and by eliminating it 
with the use of equation of motion, we get the action of the higher derivative gravity [4]. 
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It is known that the general form of the action in Eq.(l) can be transformed to a 
action in the Einstein frame as followingsfl]: We make a conformal transformation 

9„ = J^f, (2) 



and define a new scalar field as 



-1 o//'\2i 



d x = [T^)r l -KjYW- (3) 

Then, we change the scalar potential as 

U(x) = U(*)f\*) (4) 
and finally we get the action in the Einstein frame: 

S = j dix \R9[^tc R (9) + ICg^xdvX - \U{X))\- (5) 

Except for the singular case Tf^ 1 — 3(y) 2 = where we just obtain the action with no 
kinetic term of the scalar field. Hereafter, we consider the non singular case and thus we 
treat the action in Eq.(l) with T — f — 1. 

We start from the metric of the universe to be the spatially homogeneous metric as 

{dsf = a\N 2 {dt N f - a 2 {t N )[-^- 2 + r 2 (d0) 2 + r 2 sin 2 9(dcp) 2 }), (6) 

where N is the gauge function and a 2 = 2G/3n. We assume that depends only on the 
time tpf. Writting 

= (2n 2 a 2 )^(t N ), V{<j>) = 2ttVC/($), (7) 

the 1 dimensional action in Eq.(l) with / = T — 1 is expressed as follows: 

S = i J dt N [N-\-aa 2 + a 3 2 ) + N(-c?V(<P) + Ka)}, (8) 

where we have omitted the volume factor of the space. Note that the variation with 
respect to N gives us the hamiltonian constraint 

(-) 2 -0 2 -V + ^ = O. (9) 
a cr 
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An important point in our method is that the equation of motion does not depend on 
the choice of this gauge function N. In other words, we can choose this function in such 
a way that the system becomes easy to be solved. This is the main advancement from 
the method discussed in Ref. [2]. 

We are now going to change the action in Eq.(8) to the one with ordinary kinetic 
terms. All we have to do is to change variables a and <fi to z and w in such a way that 
the action is transformed to 

S = J dt N \(2N)-\-z 2 + w 2 ) - AW], (10) 

or by introducing variables z± = z ± w, we write Eq.(10) as 

S = J di JV [-(2JV)- 1 i+i_ - NW(z+, *_)]. (11) 

When we set N = 1, the cosmo logical problem turns out to be just an ordinary kinemat- 
ical problem with indefinite metric. 

In order to find new variables, we define 

a = e h{z +' z -\ T± = h±<f>. (12) 

Then the action in Eq.(8) is written as 

S = \j dt N {-^-^i + i_ + N[-e^-W{\{r + - r_)) + Ke^~\ (13) 

The change of variables which preserves the form of kinetic term is know as the two 
dimensional conformal transformation with respect to the fields r + and t_ . We therefore 
take the following choice of variables: 

t + = t + (z + ), r_ = r_(z_). (14) 

By choosing gauge function as 

N = Ne^+^d+T+d-r-, (15) 
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we find that the kinetic terms of Eq.(13) transforms to those of (11), as desired. Here d + 
and d- mean derivatives with respect to z + and Z-, respectively. To simplify Eq.(15), we 
will use the following functions in stead of r±: 



3^ 3, 



x + (z + ) = e2 T+ , X-(z_) = e2 T - (16) 
Then Eq.(15) is just written as 

N = ^Nd+x+d.x.. (17) 

In terms of x + and the original fields a and can be written as 

a=(x + x_)s, (p=-ln— . (18) 

3 £_ 

We can also find that the potential term in Eq.(lO) or (11) is given by 

W(z+,z-) = ^-d+x+d-x-ix+x-viln^] -K(x+x-)*]. (19) 

Note that the form of the potential W depends on the function x + and x_ which are 
arbitrary functions of z + and Z-, respectively, while the kinetic term takes the form of 
as seen in Eq.(ll). We utilize the freedom of choosing functions of x + and x_ to 
simplify the potential W such that differential equations are solved. 
Firstly, we consider the exponential potential with K = 0: 

[A] : V{4>) = Xe-^ (a ^ ±6) 

This case is equivalent to the Brans-Dicke theory with a constant potential and is 
known to possess a power law inflational solution a = ao(t + c) /3 (3 = as a special 
solution. The special case of a = ±2 in this case is also the special case of [H] whcih we 
discuss later. The case [H] is well studied and the analytic solution is obtained[5]. Thus, 
a = ±2 case was known to be soluble. But for other cases, the analytic solutions are not 
known before. 

Here we show that the general solution is obtained in our method. We take the ansatz 

x+ = z\ x. = z™. (20) 



When we take I and m as 

J 



'=23|- m =2Tf (21) 
with i and j are integers, then the potential W is written as 

W(z+, z.) = ^ijz^zt 1 . (22) 

36 — cr 

These are valid for a ^ ±6. When the choice of the pair is one of (1,1), (1,2), 

(2,1), (1,3), (2,2) and (3,1), the theory can be solved by the separation of variables. In 
particular, if we take i — j — 1, the potential becomes constant. With the gauge choice 
N — 1, the equations of motion are written by z + — zL — so that 

z + = a x t N + a 2 z_ = bit N + b 2 . (23) 

The Hamiltonian constraint — W = gives a restriction 

4A 

ai&i = (24) 
36 — a 2 

The scale factor a and the scalar field are obtained from Eq.(18) as functions of 
through Eq.(23). From Eq.(17) with N — 1, the time t N is expressed as a function of t 

as 

1 = oa 4 2 / + a2) ,_1 (6i*jv + h) m - l dt N) (25) 
3o — ot J 

where / and m are given in Eq.(21) with i — j — 1. In order to compare the power low 

behavior given before, we have to perform the above integration. Here we consider the 

— — 

large t N or t limit of a and 0. From Eq.(25), we have t oc t^ m ~ . Then, a = z\zl oc 
j.(i+m)/3 ^ ^--z w j 1 j c ] 1 re p ro duces the special solution stated before. 

With the ansatz in Eq.(20), some models with potentials consisting of the sum of 
exponential forms are solvable. We found the following cases with K = are solvable: 

[B] : V(<P) = Aie 2 ^ + A 2 e- 2 ^ 

We find with the ansatz in Eq.(20) 

1 3 

W = -(Xiz+ + \ 2 z-) with l = m = -. (26) 
8 4 

This case is the special case of [H] explained later so that this model was known to be 
soluble [5]. 



[C] : V((j>) = \ie-rs* + A 2 e-|^ 

We find with the ansatz in Eq.(20) 

25 5 5 

W = —{\ 1 z + + \ 2 z 2 _) with / = -, m=-. (27) 

Note that a new soluble model is derived from a model by changing <fi to -<fi. In this case, 
W is obtained by changing z + and z_ and / and m. Thus, the case 

[C] : = Aie-T * + A 2 ei* 

is also soluble and with the ansatz in Eq.(20) we find 

25 5 5 

W = — (A^_ + A 2 4) wi th Z = 7 , m = -. (28) 
144 + 4 8 

[D] : \/(0) = A ie 3 ^ + A 2 e- 3 * + A 3 

We find with the ansatz in Eq.(20) 

2 

W= -(Ai2+ + A 2 z_ + A 3 z+2-) with Z = m = l. (29) 
9 

Solutions for all of above cases can be derived following the same manner demonstrated 
for the exponential potential case. The special cases of [D], the case with Ai, A 2 > and 
A3 = — 2v / AiA 2 and the case with A 3 = are solved by Ritis et. al.[2]. 

As for the exponential potential in the case [A], a = ±6 are excluded. Here we 
consider these cases: 

[E] : V(<f>) = \e-^ 

By taking the ansatz 

x+ = e Az+ , x- = zl, (30) 

we get 

W = —A\. (31) 
18 V ; 

Thus, this case is solvable. The case with a = — 6 can be solved by exchinging x + into 

1 

X-, i.e., by taking x+ = z+ and x_ = e Az ~ . 
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Let us examine some other type of ansatz, 

x + = e Az+ , x _ = e Bz -. (32) 
There are several solvable cases under this ansatz. One is the potential 
[F] : V(<f>) = Aie^ + A 2 e^, a[3 = 36 (a ^ (3) 
In this case, with the choice 



2 + 
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B = =— |A (33) 



the potential reduces to 

W = \AB ( X, exp[2A(2 + %)z] + A 2 exp[-A(2 + . (34) 

9 V 3 a 3 / 

Another nontrivial example which cannot be solved by the separation of variable is 
the case of a Sine-Gordon potential 

[G] : V{4>) = A cos 6^0 

In this case, we take 

x + = (288)5 e -i*+, x_ = (288)5e _ 5 z -. (35) 
Then, we get the following potential: 

W = Ae _ 3*cos^iu. (36) 
In this case, we can find the following conserved currents: 

H = ~rf + ip£ + Ae-**cos^«; = 0, (37) 

1 2 1 s v^3 „z . a/3 A _z y/3 , . 

J = 2 PzPw + Q Pw + ~Y X P* e 2 sm— w + -p w e 2 cos— w, (38) 

where p z and p w are conjugate momentum of z and w, respectively. Since we found 
two conserved currents, the model is soluble. The origin of the higher dimensional con- 
served current is understood by noticing that the system can be obtained by the analytic 



continuation of the nonperiodic Toda theory [6]. Explicitly, the conserved quantities H 
and / are constructed from a Lax pair (L, M) of the nonperiodic Toda theory: That is 
H oc TrL 3 and I oc TrL 3 with the identification q± = —{z + ^=u;)/4, g 2 = an d 
?3 — ~~ (<?i + 92) where gi, q 2 and ^3 are defined in Ref.6. 

Finally, we consider the integrable models with K ^ 0. This case is in general hard 
to be solved. Except for a trivial example V = A, we can find the soluble model with the 
potential, 

[H] : V{4>) = \ x e^ + \ 2 e~^ (K ^ 0) 
Similarly to the case [B], we take the ansatz 

x+ = z^ m , x_=zl, (39) 

we obtain 

W = -(\ lZ+ + \ 2 Z- -K) (40) 

o 

so that this case is soluble. This is a famous model studeid extensively [5] and is known 
to be solble. 

It may be worthwhile to note that our result depends heavily on the fact that we 
introduced only one scalar field. In this case, we can utilize the 2d conformal transfor- 
mation with respect to a and and we found several solvable models. If we have n 
scalar fields and wish to use similar simmetry, we may need to consider n + 1 dimensional 
conformal transformation. Thus, the integrable models are quite restrictive except for 
n — 1 because the 2 dimensional conformal transformation is known to be special. Note 
also that the basic form of the soluble potential is of the exponential type. These are 
quite analogous to the integrability of 2 dimensional gravity. 

We have found some class of cosmological models can be solved in terms of the appro- 
priate choice of the gauge function and the conformal transformation. But the integrable 
models are rare and most of the physically interesting examples are not integrable. For 
example, the potential which corresponds to R 2 gravity [7] can be transformed to the 
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quartic potential by an appropriate choice of the conformal transformation. Of course, 
since we do not know the whole classification of integrable models even in the ordinary 
kinematic problem, we have an opportunity to get other integrable models. We hope 
larger class of interesting models can be found and analyzed in our method. 

Finally, we point out that for some potentials given above, the Wheeler-De Witt 
equation can be solved analytically. This equation for the case [H] ( and thus for the case 
[A] with a = ±2 and the case [B]) was solved and the analytical solution is given[5]. By 
using our transformations, we can get the analytical solutions for the cases [A], [D], [E] 
and [F] which will be given in the forthcomming paper. 
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